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The structure of the light proton and neutron rich nuclei is studied in a microscopic 
multicluster model using the stochastic variational method. This approach enables us 
to describe the weakly bound nature of these nuclei in a consistent way. Applications 
for various nuclei 6_9 Li, 7 Be, 8 B, 9 C, 9 ~ 10 Be, 9 ~ 10 B presented. The paper discusses the 
relation of this model to other models as well as the possible extension for p and sd shell 
nuclei. 

1. Introduction 

The experimentally intensively studied unstabil nuclei Jl| have challenged the theo- 
retical nuclear physicists. These nuclei exhibit a "nucleon halo", a new form of nuclear 
matter, characterized by low and nonuniform density distribution. To describe this nu- 
cleon (proton or neutron) halo care must be taken on the proper description of the nuclear 
dynamics. The halo stucture is mostly interpreted as two- or three-body phenomena, that 
is most of the theoretical models assume a stabil core and one or two valance nucleons 
@. Although one can understand various properties of the halo nuclei in these phe- 
nomenological models, it is a natural aim to go beyond these approaches. While there 
are sophisticated techniques to tackle the two- and three-body problems, the solution 
of the nuclear many-body problem is still much too complicated. The small separation 
energies a characteristic property of the halo nuclei, require extra caution. One often has 
to reproduce as small as 100 keV energy differences (for example the proton separation 
energy of 8 B or neutron separation energy of 11 Li). These tiny energies play a crucial role 
in determining the asymptotic form of the wave function and therefore the halo structure 
and the physical properties characterizing the halo nuclei (momentum distribution, inter- 
action cross section) depends on them very much. A reliable description, therefore, should 
be able to give accurate energies and correct asymptotic behaviour. Another important 



feature to be taken into account is the presence of correlations between nucleons. While 
the correlation of nucleons is a well known property in nuclear physics, it is somewhat 
surprising that there are reasonable experimental indications to assume the existence of 
correlated neutron pairs in such a low density nucleon matter as the halo of 6 He or 11 Li. 
The above metioned specific conditions (correlation, extended nonuniform density dis- 
tribution, importance of asymptotic part) make very difficult the application of models 
using single particle bases. Much of the succes of the three-body desription is due to their 
capability to treat these requirements by applying a wave function depending on relative 
variables. The antisymmetrization is, however, much more difficult in the case of relative 
coordinates and one often sacrifices the exact treatment of the Pauli principle in favor of 
simpler model, assuming structureless clusters or "core" (for example 9 Li+n+n). At this 
point, besides the validity of assuming a simplified core, another problem, the question of 
interaction between constituent particles appear. 

Our approach tries to combine the advantages of different descriptions. A microscopic 
framework is used, that is we treat nucleonic degrees of freedom thereby avoiding the need 
of knowledge of core-nucleon or nucleus-nucleus interaction. We use relative coordinates to 
have a flexible coordinate-system to describe the dynamics and to eliminate any problem 
in connection with the center-of-mass motion. Correlated basis functions is used to answer 
the challanges of the weakly bound halo structure. The spatial part of the basis functions 
are chosen to be Gaussian form to facilitate the fully analitical calculation of the marix 
elements. The Pauli principle is treated exactly. 

2. Formalism 

In our variational approach the basis functions are assumed to have the form || 



where x = (xi, Xjy_i) is a set of relative (Jacobi) coordinates, the operator A is an 
antisymmetrizer, the function 9lm l {*) represents the angular part of the wave function, 
Xsm s is t ne spin function and 4> int is the intrinsic function of the clusters. The angular part 
QlMl ( x ) is a vector coupled product of spherical harmonics of the relative coordinates. 
The spin of the clusters coupled to total spin S. 

The assumption of clusters is based on physical motivations as well as practical con- 
siderations. The clustering in light nuclei has long been known and supported by various 
experimental facts as well as numerical simulations @||. If one assumes Gaussian packet 
<^s(r) = (27/71) 3 / 4 exp{— 7(r — s) 2 } single particle states to describe a nucleus and try 
to find the optimal positions of the centers of the Gaussian packet, that is try to find 
the configurations that minimize the total energy in a realistic potential, one observes 
the formation of various clusters of nucleons (alpha, triton, 3 He, etc.). The experimen- 
tal and theoretical research devoted to unstable nuclei is also in favor of this approach 
(e.g. 6 He=a+n+n, 11 Li= 9 Li+n+n, and so on). The wave function of the clusters are 
approximated by a single harmonic oscillator shell model configuration. If more accurate 
description of the internal structure of clusters is needed one can superpose shell model 
configurations ||. Alternatively, one can divide the clusters into smaller entities, for ex- 
ample if the description of the triton with a single shell model configuration insufficient 
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one can descibe it as p + n + n. The model assumes an equal harmonic oscillator size 
parameter (3 for the clusters. 

The assumption of clusters with relatively simple description of their structure requires 
the usage of an effective nucleon-nucleon interaction, because a realistic interaction with a 
strong repulsive core would not render the clusters to be bound. The effective interaction 
used in this calculation (Minnesota interaction) contains a spin-isospin dependent central 
and a spin-orbit potential. The Coulomb interaction between protons is treated exactly. 

The matrix elements appearing in the variational equation are calculated in a fully 
analytical way ||. 

3. The stochastic variational method 

The adequate choice of the nonlinear parameters (the elements of the (N — 1) x (N — 1) 
matrix A) is very important. As the trial function contains a large number of nonlinear 
parameters and one has to superpose many basis functions to get the energy minimum a 
direct optimization of the parameters is not suitable. The number of the different spin- 
isospin and partial wave channels further complicates the choice of the basis functions. 
Moreover, the basis functions are nonorthogonal and none of them are indispensible, any 
of them can be equally represented by some other choice. This property qualifies a random 
selection of the basis function by judging on their contribution to the energy. We set up the 
basis stepwise by randomly choosing A from a preset domain of the parameter space and 
increase the basis dimension by one if the energy gain by including the randomly selected 
basis element is larger than a preset value, e. This is repeated until the energy converges. 
This random selection of the nonlinear parameters gives very accurate energy and keeps 
the size of the basis feasible. We are able to solve N = 2—7 body problems with this simple 
strategy. The results, as shown in Table 1., agree very well with those in the literature. We 
also mention that our calculations for different potentials and different physical systems 
(including Coulombic systems in atomic physics) proved to be as accurate as the most 
precise methods. The accuracy reached in these calculation justifies the application of 
this method for the field of unstable nuclei. 

The spin-isospin and partial wave channels are also randomly selected eliminating any 
bias from the construction of the wave function. Several examples show that this random 
selection gives the same percentage of different components of the wave function as a more 
direct way (e.g. choosing appropriate basis states in all possible channels and diagonalizing 
in the resulted large space) and it does not lead to a false wave function. 



4. The microscopic multicluster model 



The method has been applied for various nuclei, such as 6 8 He [jOlfTj, s 9 Li, 7 Be, 8 B 
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C, 9 - 10 Be, 9 - 10 B [0^TJ. 
Recently, we have studied the mirror nuclei 9 C and 9 Li in a microscopic a+ 3 He+p+p and 

Li mirror two- and three-body 
The calculated 



a+ 3 B.+n + n four-cluster model lOfl. The 7 Be- 7 Li and 8 B 



subsystems are also investigated with the same effective interaction [15 



ground state energies, the radii, and the densities of the nucleons are in good agreement 
with the experimental data. The magnetic and quadrupole moments, except for the 
magnetic moments of 8 B and 8 Li, are also reproduced well. The quadrupole moments of 



Table 1 

Energies and rms radii of iV-nucleon systems interacting via the Malfliet-Tjon potential. 
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9 C and 7 Be are predicted to be —5.04 e fm 2 and —6.11 e fm 2 . The microscopic multicluster 
model predicts that the neutron skin thickness is about 0.4 fm in 8 Li and 9 Li, while the 
proton skin thickness is 0.5 fm in 8 B and 9 C. Comparing to the neutron skin thickness of 
0.8 fm found in 6 He and 8 He we conclude that these nuclei do not show pronounced 
halo structure. 

The mirror nuclei 9 Be and 9 B have also been described in our model [14]. These nuclei 
are described in a three-cluster model comprising two a-particles and a single nucleon. 
The three-body dynamics of the clusters is taken into account by including both of the 
possible arrangements, (aa)N and (Na)a, and by using all the relevant partial waves of 
the relative motion of the clusters. The ground state of 9 Be, the only particle-bound state 
in the spectra of these nuclei, is calculated by using the stochastic variational method, 
while the other particle-unbound states are studied by the complex scaling method. The 
calculated spectra of 9 Be and 9 B are compared with experiment in Table 2. The theoretical 
level sequence in 9 Be has a good correspondence with the observed spectrum. The second 

3 
2 



resonance is obtained at 4.3 MeV excitation energy. The other calculations |T6Hl8 
also predict the § state. Although no such state is cited in Ref 



19 



the calculated 

resonance may correspond to the state at 5.59 MeV mentioned in Ref. [20|. We get 
two broad overlapping resonances with | 
the conclusion of the recent experiments 



and | 



at about 6.5 MeV. This agrees with 



20,21 



We could not find a resonance with | 



around 8 MeV excitation energy in accordance with Refs. 
is parenthetically quoted in Ref. 



20,21 



19 



Instead of this a | 



although such a state 
resonance is obtained at 7.9 

MeV, which agrees with the result of Refs. [IB, 171. The spectrum of 9 B is less known 
experimentally compared to that of 9 Be. The calculated spectrum is similar to the one of 
9 Be. We can predict the energy and the width of several resonances in 9 B with the same 
accuracy as the case of 9 Be. For example, our calculation predicts a missing ~ state of a 
1 MeV width at the excitation energy of 2.43 MeV, which is in agreement with the result 
of a recent 9 Be(p, n) reaction [^2] that located the \ state at 2.83 MeV. Although no 
definitive spin assignment is made to the state at 2.788 MeV excitation energy [jn|, our 
calculation supports gnment rather than | . The first excited 1/2+ state was 

not localized in the present study. 



Table 2 

Energies and widths of the unbound states in 9 Be and 9 B. The energy is from the three- 
body threshold. The spin and parity of the 3.065 MeV state of 9 B is assumed to be 
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The electromagnetic moments and the rms radii of proton, neutron, and nucleon, assum- 
ing pointlike nucleons, are included in Table 3. Bare operators are used in the calculation. 
The theory is found to reproduce the data very well. The fact that the total reaction cross 
section is reproduced well also supports that the calculated ground state density is re- 
liable. The 1/2+ -> 3/2 _ El transition and the 5/2~ -> 3/2 _ E2/M1 transitions are 
studied by treating the excited states as quasibound states. The calculated transition 
rates, shown in Table 3., are in good agreement with the experiments. 



Table 3 

Radii and electromagnetic properties of 9 Be. The reduced matrix elements are given in 
Weisskopf units. The bare-nucleon charges and ^-factors are used in the present calcula- 



tion. The effective charges were used in the shell model calculation of Refs. pl| and [23 



to calculate the quadrupole moment and the E2 strength. See text for the B(E1) value 
of the present calculation. 
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The fact that the present calculation reproduced all the data very well strongly supports 
that the three-cluster model is quite appropriate for describing the structure of 9 Be and 
9 B, provided that the three-body dynamics is treated properly in the calculation. Because 
the ground state and the | ,2.43 MeV state are described well by the present model, the 
(3 decay of the 9 Li ground state to these states is expected to further test the accuracy 
of their wave functions or an available wave function of 9 Li. The experimental value of 
log/t for the (3 decay to the 9 Be ground state is about 5.31 [|19|,^4[], indicating that the 
/9-decay matrix element is fairly suppressed despite the allowed transition. The weak (3 
decay is ascribed to the fact that the spatial symmetry of the main component of 9 Be is 
different from that of 9 Li [26]. In fact the Gamow- Teller (GT) matrix element to any state 



of 9 Be, if it is described by the a + a + n three-cluster model, always vanishes regardless 
of the wave function of 9 Li. To explain the weak f3 decay we have to admix a small 
component which is not taken into account in the three-cluster model. This possibility 
has been investigated by including the distortion of the a-particle into t + p and h + n 
configurations. The calculated log/t value turns out to be 5.60, indicating that we are on 
the right track. By being able to accommodate such distortion into the model consistently 
we have exemplified a unique advantage of the microscopic multicluster model. 

We have investigated the question of whether or not some high isospin excited states of 
stable nuclei have extended halo-like structure. For this we have focused on the + state 
of 6 Li at 3.563 MeV excitation energy, which is the T 2 =0 member of the isospin triplet 
together with the ground states of 6 He and 6 Be. This indicates that the spatial structure 
of 6 Li(0 + ) is very similar to the 6 He ground state which is known to have two-neutron 



halo structure. 

We have done microscopic calculations for 6 Li(0 + ) with a+p+n three-body model by 
allowing the isospin mixing due to the Coulomb potential. A large number of channels was 
taken into consideration in the calculation. The isospin mixing was found to be moderate 
and does not destroy the isobaric analogue concept proposed in ref. |25| even near the 



three-particle threshold. The accuracy of the wave function was tested by calculating the 
electron scattering form factors. The calculated form factors are in fair agreement with 
experiment. 

The neutron and proton density distributions of 6 Li(0 + ) are similar and more extended 
up to larger distance than the 6 He density. The matter size of 6 Li(0 + ) is calculated to 
be 2.73 fm, which is by 0.2 fm larger than the 6 He size. Our analysis strongly supports 
that the 3.563 MeV + state of 6 Li has spatially extended halo-like structure formed by 
the neutron and proton outside the a-particle. The inelastic proton [^J and pion |2"8[ 
scatterings to this state show strong anomaly in the sense that any theoretical analysis 
leads to outstanding discrepancy from experiment. We think the consideration of the 
halo-like structure in this state is needed in such an analysis. 

Further development of radioactive beam facilities will open up a spectroscopic study 
on high isospin excited states of stable nuclei as well as unstable nuclei near the drip-line 
and give us valuable information on the isospin impurity and the generalization of the 
neutron-halo concept. 



5. Comparision to other models 

As we have shown in a number of examples, the microscopic multicluster model pro- 
vides us with a good and consistent description of light nuclei including unstable nuclei. 
The reason for success is that it can duly take into account the dynamical correlation be- 
tween the nucleons as well as the asymptotic behavior characteristic to the weakly bound 
system. To treat the relative motion of the clusters flexibly enough, the trial function for 
the relative motion is chosen to be correlated Gaussians . The complexity of the wave 
functions is increased further if different cluster partitions are included in the model space. 

The direct comparison of our wave function and model space to that of shell model 
or even to other versions of cluster models (which use different trial functions) is very 
complicated. Some way for comparision, hovewer, would be inevitably necessary to have 
deeper understanding of applicability and limitations of different models. It would be 
useful to have a simple measure to compare various types of wave functions. We have 
recently proposed to use the occupation probability of the number of total oscillator 
quanta as a possible candidate as a tool of comparison p3f . The probability Pq of a 
definite number of total HO quanta Q can be obtained by calculating the expectation 
value of the operator O 

= ±- [ 27T d9 exp(z0 ( Y.P l [AhoCO " |] " Q)\ (2) 

JO ■ £ 

where H-ro(i) is the 3-dimensional HO Hamiltonian divided by Hoj = —7. 



Pi projects out either proton or neutron. It is set the unit operator when one calculates 
the number of total quanta occupied by both protons and neutrons. The advantege of 
this formalism is that the evaluation of the matrix element of this operator is very simple 
in cluster models and it is trivial in shell model. Table 4. shows the examples for a pair of 
the mirror nuclei 7 Li- 7 Be, 8 Li- 8 B, and 9 Li- 9 C. The result for 6 He and 9 Be is also included 
in the table. The value of 7 is set 0.17 fm -2 (tkv =14.4 MeV). The probabilities are given 
as a function of Q exc = Q — Q m i n , where Q m i n is the minimum number of HO quanta for 
the lowest Pauli-allowed configuration. The lowest Ohu component is around 50-60 % for 
most cases and the sum of 0, 2, and 4hu components accumulates to about 90 %. The 
admixtures of higher components than Q exc = 4 are significant in the ground states of 8 B 
and 9 Be and also in the ground state of 6 He. 

6. Extension to larger systems 

The application of the SVM on correlated Gaussian basis for more than N=6-7 body 
system is difficult as both the partial wave expansion and the calculation of the matrix 
elements of the Hamiltonian become too complicated and computer time consuming. Let 
us try to use some other type of basis functions to avoid these problems. The Gaussian 
packet (p2(r) = (2z//7r) 3 / 4 exp{— 7(r — s) 2 } functions are also often used in few-body and 
few-cluster calculations. These functions do not belong to a particular orbital angular 
momentum and if one uses them as single particle basis functions in a Slater determinant, 
the calculation of the matrix elements will not cause serious problems. In a variational 
calculation the basis functions take the form: 



The wave function of the system is approximated by linear combinations of these Slater 
determinants. States of good angular momenta can be obtained by letting an angular 
momentum projection operator act on the full wave function. 

The adequate choice of the nonlinear parameters {i>i, Sj}^ of the Slater determinants is 
very important. These nonlinear parameters can be selected by the stochastic variational 
method (SVM) or by an appropriate direct optimatization. To keep the cost of the 
optimatization low only a set of the nonlinear parameters in a given Slater determinant 
was optimized at a time, while those of the others were kept fixed. 

The energy quite slowly converge on this uncorrelated basis and for few-body systems it 
fails to reach the same value as on the angular momentum projected Gaussian (see Table 
5.). For heavier systems, where the angular momentum projected Gaussian basis is not 
feasible, the results on shifted Gaussian basis is quite close to those of other methods such 
as the Integrodifferential Equation Approach (IDEA) or the Hartree-Fock-Bogoljubov 
(HFB) method. 

The Antisymmetrized Molecular Dynamics (AMD) || and Fermionic Molecular Dy- 
namics (FMD) || methods use the shifted Gaussian basis functions and often approxi- 
mate the wave function of a multinucleon system by only one Slater-determinant (K = 1 
in Table 5.). Table 5. shows that this approximation can be insufficient, and the linear 
combination of Slater determinants considerably lowers the ground state energy. 
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Table 4 

The occupation probability of the number of harmonic-oscillator quanta for microscopic 
multicluster-model wave functions. The probabilities for nucleons, protons, and neutrons 
are given in % in the upper, middle, and lower rows, respectively, as a function of oscillator 
excitations. Asterisk indicates the probability of less than 1 % and dashed line represents 
vanishing probability. The average number of oscillator excitations is given in the column 
labeled < Q cxc >. The details of the wave functions are referred to Ref. [12] for 6 He , to 
Ref. [13] for 7 Li, 8 Li, 9 Li, and 9 C, and to Ref. [14] for 9 Be. 
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Table 5 

The binding energies (in MeV) of different A-nucleon systems interacting via the Volkov 
potential (m = 0.6) 
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-8.46 


4 He 


-27.92 


-29.75 


SVM 


-30.42 


6 He 


-23.59 


-29.38 


SVM 


-31.82 


8 Be 


-52.60 


-57.09 






16Q 


-1100.1 




IDEA 


-1101. 



7. Summary 

In summary, we presented the new results of the microscopic multicluster model using 
the stochastic variational method. The stochastic variational method provides us with 
accurate energies and wave functions. We have determined physical quantities of interest 
(proton, neutron and matter distributions and radii, magnetic and quadrupole moments, 
electromagnetic transition rates, beta-decay probabilities, momentum distribution of frag- 
ments, interaction cross sections, spectroscopic amplitudes) for various light nuclei. We 
have discussed the relation of our approach to other models. Further applications for light 
nuclei, for example 10 Be, 11 Be, n Li are under way. 
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